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^ ■ Abstract 
D 

["t I ■ Multi-rate asynchronous sub-Nyquist sampling (MASS) is proposed for wideband spectrum sensing. 

' Corresponding spectral recovery conditions are derived and the probability of successful recovery is given. 

Compared to previous approaches, MASS offers lower sampling rate, and is an attractive approach for 
cognitive radio networks. 
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OO 

■ I. Introduction 

O 

' The radio frequency (RF) spectrum is a limited natural resource, which is currently regulated by 
government agencies. The primary user (PU) of a particular spectral band has the exclusive right to use 



^ ! that band. Nowadays, on the one hand, the demands for the RF spectrum are constantly increasing due to 

the growth of wireless applications, but on the other hand, it has been reported that the spectrum utiUzation 
efficiency is extremely low. Cognitive radio (CR) is one of the promising solutions for addressing this 
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spectral under-utilization problem |IT]. An essential requirement of CRs is that they must rapidly fill in 
spectrum holes (i.e., portions of the licensed but unused spectrum) without causing harmful interference 
to PUs. This task is enabled by spectrum sensing, which is defined as a technique for achieving awareness 
about the spectral opportunities and existence of PUs in a given geographical area ||2]-||4]. 

CR with a "broader" spectral awareness could potentially exploit more spectral opportunities and 
achieve greater capacity ||5]. Wideband spectrum sensing techniques therefore have attracted much at- 
tention in the research on CR networks. In |)6l, a wavelet-based approach was studied for performing 
wideband spectrum sensing. It provides an advantage of flexibility in adapting to a dynamic wideband 
spectrum. Quan et al. Q, HI proposed a multiband joint detection (MJD) approach for detecting the 
primary signal over multiple frequency bands. It has been shown that MJD performs well under practical 
conditions. In |9 |, a filter-bank method was presented for sensing the wideband spectrum in a multicarrier 
communication system. The filter-bank system has been shown to have a higher spectral dynamic range 
than other traditional power spectrum estimation methods. Furthermore, compressive sensing (CS) lITOl . 
inn technologies were introduced to implement wideband spectrum sensing in |[T2l - |[T7l . Notably, these 
techniques take advantage of using sub-Nyquist sampling rates for signal acquisition, instead of the 
Nyquist rate, leading to reduced computational burden and memory requirements. 

In this paper, we develop a multi-rate asynchronous sub-Nyquist sampling (MASS) system to perform 
wideband spectrum sensing. As the spectral occupancy is low, spectral aliasing (generated by sub- 
Nyquist sampling) is induced in each sampling branch to wrap the sparse spectrum occupancy map 
onto itself. The sensing requirements are therefore significantly reduced, i.e., the proposed MASS system 
has superior compression capability compared with the Nyquist sampling system. Considering the sub- 
Nyquist sampling in MASS, we then determine the recovery conditions under which the full wideband 
spectrum can be successfully reconstructed by using CS analysis. Compared with existing wideband 
spectrum sensing approaches, MASS has lower implementation complexity, higher energy efficiency, 
better data compression capability, and is more applicable to CR networks. 

The rest of the paper is organized as follows. We briefly introduce the CS-based sensing scheme 
in Section lUl We then propose the MASS system in Section Hill Simulation results are presented in 
Section Jvl followed by conclusions in Section Ivl 
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Fig. 1. Diagram of CS-based spectrum sensing when using the spectral domain energy detection approach. 

II. Problem Statement 

In this paper, we assume that all CRs keep quiet during the spectrum sensing interval as enforced 
by protocols, e.g., at the medium access control (MAC) layer |7|. Therefore, the observed signal at 
a CR arises only from PUs and background noise. Suppose that the continuous-time signal Xc{t) is 
received at a CR, and the frequency range of Xc{t) is ~ (Hz). If the signal Xc{t) is sampled at the 
sampling rate /s (Hz) for an observation time T, the sampled signal can be denoted by x[n] = Xc{n/ fs), 
71 = 0, 1, • • • , A^— 1, (in vector form x £ C^^^), where = fsT is assumed to be an integer. The discrete 
Fourier transform (DFT) spectrum of x can be calculated hy X = Fx, where F denotes an A^-by-A^ DFT 
matrix F[/i, n] = e-i27rn/i/7V ^vhich j = The difficulty arises when we consider the Shannon- 

Nyquist sampling theorem, which requires the sampling rate to be at least twice the bandwidth of the 
signal, i.e., /s > 2W. On the one hand, we would like to realize "broader" spectral awareness at CRs (i.e., 
larger W), but on the other hand, the higher sampling rate will result in excessive memory requirements 
and prohibitive energy costs. This dilemma has motivated researchers to look for technologies to reduce 
the sampling rate /s while retaining W by using CS theory ifTOl . lITTI . 

CS theory indicates that if the signal is sparse in a suitable basis, it can be exactly recovered from 
partial measurements. By "sparse" we mean that the signal can be represented using a few coefficients in 
some basis. Based on the fact of spectral sparseness [13], it is reasonable to sample the signal Xc{t) at a 
sub-Nyquist rate, reconstruct the spectrum X from partial measurements, and perform spectrum sensing 
using the reconstructed spectrum X. A typical spectrum sensing approach is spectral domain energy 
detection ifTSl . As depicted in Fig. [TJ this approach extracts the reconstructed spectrum in the frequency 
range of interest, e.g.. A/, and then calculates the signal energy in the spectral domain. The output energy 
is compared with a detection threshold (denoted by A) to decide whether the corresponding frequency 
band is occupied or not, i.e., choosing between hypotheses Tii (presence of PUs) and Tio (absence of 
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Fig. 2. Schematic illustration of the multi-rate asynchronous sub-Nyquist sampling system in one CR node. The wideband 
filter has bandwidth W . 

PUs). It is clear that the performance of sub-Nyquist-based spectrum sensing will highly depend on the 
recovery quality of the spectrum. In this paper, we will present a novel system, i.e., MASS, to sample the 
signal using sub-Nyquist sampling techniques, while enabling the spectrum X successfully recovered. 



III. Multi-rate Asynchronous Sub-Nyquist Sampling 
A. System and Signal Model 

Suppose that a CR has v sub-Nyquist sampling branches as shown in Fig. |2] The wideband filter prior 
to the samplers removes frequencies outside the spectrum of interest, and is set to have bandwidth W . At 
the i-th branch, the low-rate sampler samples the received signal at the sub-Nyquist rate fi (Hz). In the 
observation time T (second), the numbers of samples in these v sampling branches are Mi, M2, • • • , My, 
respectively, where Mj = fiT (Vi G [1, i;]). In addition. Mi, M2, • • • , M^ are chosen to be different prime 
numbers that are of the order of ^/N, i.e., Mj ~ 0{-s/N), by controlling the sampling rate /j (Vi € [1, v])- 
The DFT spectrum of the sampled signal is then computed by applying the fast Fourier transform (FFT) 
to the samples in each branch. After that, these DFT spectra are used to reconstruct the wideband 
spectrum. We will employ an energy detection approach, and so we are interested in reconstructing only 
the magnitude of the spectrum, i.e., \X\, resulting in the spectral magnitude estimate \X\. 

Suppose that the received signal Xc{t) is of finite support and absolutely summable. Using the sub- 
Nyquist rate /j < 2W, we obtain the sampled signal yi[m\ = Xc{m/ fi) = Xc(^), m = 0, • • • , Mj — 1. 
The DFT spectrum of yi is then calculated by Yi = F^yi, where denotes the Mj-by-Mj DFT matrix. 
The DFT spectrum of iji is related to the continuous-time Fourier transform of Xc{t) by |[T9l 



(1) 



Z=— 00 
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where Xc{f) = Xc(t)e~-'^'^*^ dt is the Fourier transform of Xc{t). In other words, the DFT spectrum 
of the sampled data is a sum of the shifted Fourier transform spectrum of the continuous-time data. 

Furthermore, if Xc{t) is sampled at or above the Nyquist rate, i.e. fs = Y — sampled signal 

can be written as x[n] = Xc(n//s) = Xc(^), n = 0, 1, • • • , — 1. The spectrum of x will be related to 
the continuous-time Fourier transform of Xc{t) by X{f ) = /s X]£-oo "'^c(/ + ^/s)- As the signal Xc{t) 
is band-limited to W and the sampling rate /s > 2W, there will be no spectral aliasing phenomena in 
X{f); thus, we can rewrite this relationship by X{f) = /sXc(/), V/ G [—^, ^] ■ Because X{f) has all 
information in [— we assume that X{f) is zero everywhere except / € [— Substituting 
X{f) = fsXcif) into (dl), we can obtain 



Is , 

i=—oo 



w w 



(2) 



which has the discrete form of 
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(4) 



where [a\ is the floor function that gives the largest integer not greater than a, and 5[n\ denotes the Dirac 
delta function that is zero everywhere except at the origin, where it is one. In matrix form, we write 

M, 



N 



^iX 



(5) 



where the elements of *i G j^m.xAt (^j^. ^ jy-^ represented by *i [m + |_^J + 1, n + |_yJ + l] = 

X];^_oo S[n — {m + /Mj)] . It can be seen that in each column of there is only one non-zero element 
and this has the value of one. In each row of the number of non-zero elements is at most \^~\ (the 
ceiling function gives the smallest integer not less than -^), which is the undersampling factor. 



B. Multi-rate Sub-Nyquist Sampling 

Due to the sub-Nyquist sampling in each sampling branch, we should consider the effects of spectral 
aliasing. However, when the spectral sparsity k <^ N and the sampling rate satisfies Mi ~ 0{y/N), the 
probability of signal overlap is very small. The reader is referred to Appendix A for the proof. In such 
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a scenario, we concentrate on considering two cases: no signal at a particular value m and one signal at 
m. If there is no signal overlap, the following equation holds by using ([S]): 



\YA 



N 



Mi 
N 



-^i\X\ 



(6) 



where the last equation holds because the elements of $j are either zeros or ones, while each frequency 
bin of Yi has no signal overlap from X. 

Furthermore, since all samplers observe the same spectral magnitude, i.e. \X\, a concatenated equation 
relating \X\ to the partial measurements can be formed by using 1^, i.e.. 



JL\v 



\X\ = 



(7) 



where $2, • • • , *t, are disjoint submatrices of $ G M(^i=i A^>)xAr 

Remark 1: The proposed system does not require exact synchronization between sub-Nyquist samplers. 
We merely require that the time offsets between sampling branches are sufficiently small so that the 
observed spectral magnitudes are quasi-stationary. This is because, after the signal sampling and the 
spectral recovery, we employ the energy detection approach for performing spectrum sensing, and the 
use of energy detection only requires the spectral magnitude information. Due to this advantage, MASS 
can be applied to cooperative CR networks in which each CR requires one sub-Nyquist sampler and 
the spectral environments of different CRs are similar. In such scenarios, cooperative CRs are required 
to send the compressed magnitude information |yi|, • • • ,11^1 to a fusion center for spectral recovery by 
using ©. 

To reconstruct the spectral magnitude (i.e., \X\) using (|7]|, we should carefully choose the measurement 
matrix We employ the mutual coherence fTOl to evaluate the suitability of 

Definition 1: Let $ be expressed as $ = [(pi </>2 • • • '/'a'], where (p^ denotes the h-th column of the 
matrix Then the mutual coherence of the matrix $ is given by 



u = max < 01,1 

l^zGll,N] 



> 



(8) 
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where (pi = i' denotes the £2 normaUzed column. 

\\<l>ih 

The aim is to keep /i to a minimum that allows the spectrum to be successfully recovered. When 
the conditions in the following proposition are satisfied, the mutual coherence of $ in MASS can be 
determined by the number of sampling branches v. 

Proposition 1: Let v samplers observe the spectral magnitude in the same observation time T and 
generate v measurement vectors, i.e., |yi|, jl^l; ■ " " A'^vl- If the lengths of the measurement vectors are 
different primes, Mi,M2,--- ,M^, that satisfy M/M^ > N, \/l,z G [l,v],l ^ z, then the mutual 
coherence of the measurement matrix $ is given by 

/i = max I < (^/, (/)^ > I = -. (9) 

The proof of Proposition 1 is given in Appendix B. 

Remark 2: Donoho and Elad |[20l have proven that when the mutual coherence of the measurement 
matrix satisfies fi < 2k-i ^ A:-sparse signal can be successfully recovered. Thus, we know that if 
/i = i < i-e-> V > 2k — 1, the spectral magnitude \X\ can be exactly reconstructed. Proposition 1 

illustrates that the number of samples in each sampler (or CR) is of the order of ^/N. When Mi ~ 0{VN), 
at least 2k sampling branches are required to recover the spectral magnitude. To summarize, MASS needs 
the total number of measurements to be X]i"=i ~ 0{ky/N). 

Proposition 1 requires that the number of sampling branches w > 2A; — 1 in order to reconstruct the 
spectrum. Actually, if the number of sampling branches v is less than 2k, we could still reconstruct 
the spectrum with high probability. Next, we will investigate the probability of successful spectral 
reconstruction when we have an insufficient number of sampling branches, i.e., v < 2k. 

Proposition 2: Let v samplers observe the spectral magnitude in the same observation time T and 
generate v measurement vectors, i.e., jl^l, • • • ,\'^v\- If the lengths of the measurement vectors are 
different primes and satisfy MiM^ > N, z G [1, u], / 7^ 2, then the probability of successful spectral 
reconstruction is at least 1 — J2^=i jf- 

The proof of Proposition 2 is given in Appendix C. 

Remark 3: Proposition 2 shows that the probability of successful spectral reconstruction increases as 
either the number of sampling branches or the sub-Nyquist sampling rates increase. When v = 2k — l, the 
probability of successful spectral reconstruction is at least 1 — XliLi W'- Further, as Mi ~ 0{VN) where 
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N = /sT > 2WT, we know that, with "broader" bandwidth W, we could achieve higher probability of 
successful spectral reconstruction. 



In simulations, we consider that there are v sub-Nyquist sampling branches, each of which is equipped 
with a single low -rate sampler. In the i-th sampling branch, the received signal Xc{t), as defined below, 
is sampled at a sub-Nyquist rate /j over an observation time T. 



where sinc(2;) = ^^'^^^^ , A denotes a random time offset between sampling branches, and z{t) is unit 
additive white Gaussian noise (AWGN). For the fading case, we generate v different received powers Ei 
according to the fading environment, and regenerate them for the next observation time. The received 
powers El do not change within T, but vary randomly from branch to branch subject to the channel 
fading. The wideband signal Xc{t) consists of Nf, = 30 non-overlapping subbands, whose bandwidths 
Bi = 0.5 ~ 5 MHz, with carrier frequencies /; = ~ 20 GHz. Since the signal has a bandwidth of 
W = 20 GHz, if it were sampled at the Nyquist rate for T = 2 /is, the number of Nyquist samples 
would be = 2WT = 80, 000. However, in MASS we use v = 22 sampling branches to sample the 
wideband signal with different sub-Nyquist rates, where Mj ~ 0{^/N). Specifically, we select the first 
prime Mi ^ a^/N (a > 1) and its u — 1 neighboring and consecutive primes. The spectral magnitude 
|X| is reconstructed by applying the compressive sampling matching pursuit (CoSaMP) algorithm [21 J 
to ([7]). The spectral occupancy status is then decided by using the energy detection algorithm on the 
reconstructed spectrum. 

Fig. [3ta) depicts the influence of the sparsity level k and the compression ratio ^ (M = ^'=J ) on 
the detection performance of MASS. It illustrates that the lower the sparsity level, the better the detection 
performance. In addition, a higher compression ratio will lead to a lower probability of false alarm and a 
higher probability of detection. In Fig. Sb), we can see that the proposed system enables us to perform 
wideband spectrum sensing using fewer measurements (thus lower sampling rate). This is because, for 
the compression ratio below 0.2, the proposed system can achieve smaller reconstruction mean square 
error (MSE) when compared with both the CS-based approach and the finite rate of innovation approach. 

Febraary 8, 2013 DRAFT 



IV. Simulation Results 
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Fig. 3. (a) Influence of the sparsity level and the compression ratio on the detection performance of MASS, when the signals 
from PUs suffer AWGN channels with average SNR= 10 dB. (b) Comparison between the proposed system and the existing 
approaches when the compression ratio varies. 




Number of sampling branches Probability of false alarm 
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Fig. 4. (a). Comparison of spectrum recovery performance for synchronous samplers and asynchronous samplers, with A'^;, = 30 
and Af, = 50 over AWGN channels with average SNR= 10 dB. (b). Performance of MASS over AWGN, Rayleigh, and shadow 
fading channels, with the number of sub-bands Ni, = 30, the compression ratio of 49.36%, and the standard deviation of shadow 
fading 5 dB. 

Fig. IHa) shows the effect of imperfect synchronization between sampUng branches. Compared with a 
reference clock, the asynchronous sampUng branches have time offsets in the range of ~ 0.8 /US, 
while the total observation time is 2 ^us. It is evident that the detection performance of the asynchronous 
samplers is roughly the same as that of the synchronous samplers. Fig. Ua) also illustrates that with 
more sampling branches, better spectrum sensing performance can be achieved. This is because with 
more sampling branches, a higher probability of successful spectrum recovery can be obtained. 

In view of Remark 1 and the MASS'S robustness against the imperfect synchronization in Fig. @ta). 
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the proposed system could also be applied to CR networks in which each CR only has one sampling 
branch. In such a scenario, different branches could face independent and identically distributed (i.i.d.) 
fading channels. In Fig. Ub), we analyze the detection performance of MASS over AWGN, Rayleigh, 
and shadow fading channels. When the signal-to-noise ratio (SNR) is zero decibels (dB), the performance 
of MASS over fading channels is roughly the same as that over non-fading AWGN channels. This is 
because the strength of the signal is mostly masked by the noise. By contrast, the detection performance 
of MASS over AWGN channels exceeds that over fading channels when SNR= 5 dB. In addition, it is 
found that the performance of MASS over shadow fading channels is the poorest, in comparison to the 
case of AWGN and Rayleigh fading channels. Nonetheless, even over shadow fading channels, MASS 
has a probability of nearly 80% of detecting the presence of PUs when the probability of false alarm is 
10%, with the compression ratio ^ = 49.36%. On the contrary, Nyquist sampling systems, e.g., wavelet 
detection and multiband joint detection, must use at least N measurements. 



TABLE I 

Comparisons between wideband spectrum sensing techniques. 



Technique 


Compression 
Capability 


Sampling 
Scheme 


Implementation 
Complexity 


Wavelet detection 


no 


Nyquist 


low 


Multiband joint detection 


no 


Nyquist 


high 


Filter-bank detection 


no 


Nyquist 


high 


CS-based detection 


yes 


sub-Nyquist 


medium 


Proposed MASS 


yes 


sub-Nyquist 


low 



In Table IH we can see that, compared with Nyquist sensing techniques, MASS has outstanding 
compression capability due to multi-rate sub-Nyquist sampling. A similar- system that has compression 
capability is the CS-based system. However, the state-of-the-art CS-based approaches require pseudo- 
random sequence generators as compression devices. In order to exploit spatial diversity in CR networks, 
each CR node should have a separate compression device, and transmit both compressed data and its 
measurement matrix to the fusion center. Even if measurement matrices are pre-stored, synchronization 
issues should still be considered because the un-synchronized measurement matrices could lead to false 
spectral recovery. In contrast, no separate device is required for generating the measurement matrix in 
MASS while the synchronization requirements are relaxed as shown in Fig. Ufa). Thus, compared with 
the CS-based approaches, MASS has lower implementation complexity. 
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V. Conclusions 

In this paper, we have proposed a novel wideband spectrum sensing system, i.e., MASS. In this scheme, 
parallel low-rate samplers are used to sample the wideband signal at different sub-Nyquist rates. We 
have derived conditions for recovering the full spectrum by using CS analysis. The successful recovery 
probability has also been given. Simulation results have shown that MASS has superior compression 
capability compared with Nyquist sampling systems. Unlike other sub-Nyquist sampling systems, e.g., 
CS-based systems, MASS has been seen to be robust against lack of time synchronization and to have 
excellent performance in fading/shadowing scenarios. In summary, the proposed MASS not only has 
the benefits of low sampling rate, high energy-efficiency, and compression capability, but also is more 
amenable to implementation in CR networks in the presence of fading/shadowing. 

Appendix A 

Probability of Signal Overlap While Sub-Nyquist Sampling 

As the Nyquist spectrum is A:-sparse, the probability of bin n belonging to the spectral support is 
P = jj. On letting q denote the number of spectral components overlapped on bin m, the probability of 
no signal overlap can be given by using ©-(IS]! as 

Pr(g < 2) = Pr(g = 0)+Pr(g = 1) = (l-P)^^V f P(l-P)^^^~\ (11) 



1 



Substituting P = jr into (fTTI) while using Mi = \fN , we obtain 



Pr(o < 2) = ( )".H ( ) ^'^ = ^ ^ -. (12) 

^ ^ N ^ Mi N ' N-k 

It can be tested that Pr(g < 2) approaches 1 when k <^ N. Thus the probability of signal overlap 
approaches zero. □ 



Appendix B 
Proof of Proposition 1 

Let and (ph denote the n-th and the h-th (h 7^ n) normalized column of the measurement matrix 
respectively. Using (|5), they will have non-zero elements on the row indices mi = \n\ (^j^) and 
9i = 1^1 mod (Ml) ill the l-th sub-matrix, i.e., Obviously, mi could equal gi when n| is a multiple of 
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Ml. To avoid this happening again in another sub-matrix, e.g., $2, we could let Mi and be different 
primes such that MiMz > N. This is because, if Mi and Mz are different primes and MiMz > N, 
\h — n\ cannot be a multiple of both Mi and Mz- Thus, the maximum correlation of different columns 
in $ exists only when mi = gi in a single sub-matrix ^i, in which case the mutual coherence /i = i. □ 

Appendix C 
Proof of Proposition 2 



Based on Doob's maximal inequality II221I . the following inequality holds: 



Pr(^ > x) = Pr ( max | <?!<2 > | > x I < 



E 



< (Ph<Pz > 



(13) 



Let us define the event 17/ such that mi equals gi as defined in the Appendix B, and D/ = to 
be the undersampling factor in the l-th sampling branch. The probability of the event $7/ occurring is 
Pr(il/) = '^^^f °' = ^'^^^'j^^'j^"^-* = < jj-^, where denotes the binomial coefficient. Then the 

expected value of < (pj , (ph > becomes 



E 



< (Pl,(Pz > 



V ^-^ V ^-^ Mi 



Replacing x in ([T3] ) by x = 2k-\ ' following inequality can be obtained by using (fT3l)-(fT4l): 

Pr > ^) < i2k - 1)E [< 0,, 0. >]<'-^±^. (15) 



Using the above inequality, we obtain 

Considering the successful recovery condition in [,20 J . i.e., /i < 2fcrT> we complete the proof. □ 
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